Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN WF$L style file vl.4) 



5-1 



Or 

6 

H 



CO 

c3 



Comparing the SFI Peculiar Velocities with the PSC2; 
Gravity Field: a VELMOD analysis 

E. Branching 2 , W. Freudling 3 ' 4 , L.N. Da Costa 4 , C.S. Frenk 5 

R. Giovanelli 6 , M.P. Haynes 6 , J.J. Salzer 7 , G. Wegner 8 , I. Zehavi 9 



1 Kapteyn Institute, University of Groningen, Landleven 12, P.O. Box 800, 9700, Groningen, The Netherlands. 

O ' 



2 Dipartimento di Fisica dell'Universitd degli Studi "Roma THE", Via della Vasca Navale 84, 1-00146, Roma, Italy. 

3 Space Telescope - European Coordination Facility, Karl-Schwarschild Strasse 2, 85748, Garching, Germany. 
CIh, 4 European Southern Observatory, Karl-Schwarschild Strasse 2, 85748, Garching, Germany. 

' 5 Department of Physics, University of Durham, South Rd., DH1 3LE, Durham, U.K. 

\J^} • 6 Center for Radiophysics and Space Research and National Astronomy and Ionosphere Center. 

CN | Cornell University, Ithaca, NY 14853, USA. 

7 Astronomy Department, Wesleyan University, Middletown, CT 06457, USA. 
■ 8 Department of Physics and Astronomy, Darmouth College, Hanover, NH 03755, USA. 

9 NASA/Fermilab Astrophysiscs Group, Fermi National Accelerator Laboratory, Batavia, IL 60510-0500, USA. 

m ' 
o 



1 February 2008 



ABSTRACT 



We compare the peculiar velocities derived from the /-band Tully-Fisher 
(TF) relation for 989 field spiral galaxies in the SFI catalogue with the pre- 
dicted velocity field derived from the IRAS PSCz galaxy redshift survey. We 
assume linear gravitational instability theory and apply the maximum like- 
lihood technique, VELMOD (Willick et al. 1997b), to SFI galaxies within 
cz< 6000 kms -1 . The resulting calibration of the TF relation is consistent 
with a previous, independent calibration for a similar sample of spirals re- 
siding in clusters. Our analysis provides an accurate estimate of the quantity 
Pi = ft^f/bj, where bj is the linear biasing parameter for IRAS galaxies. Us- 
ing the forward TF relation and smoothing the predicted velocity field with 
a Gaussian filter of radius 300 kms -1 , we obtain /3j = 0.42 ± 0.04 (l-cr un- 
certainty). This value, as well as other parameters in the fit, are robust to 
varying the smoothing radius to 500 kms -1 and splitting the sample into 
spherical shells in redshift space. The one exception is the small-scale velocity 
dispersion, a v , which varies from ~ 200 kms -1 within cz^q — 4000 kms -1 to 
~ 500 kms -1 at larger distance. For /3j ~ 0.42, the residuals between the TF 
data and the PSCz gravity field are uncorrelated, indicating that the model 
provides a good fit to the data. More generally, a \ 2 statistic indicates that 
the PSCz model velocity field provides an acceptable (3a) fit to the data for 
0.3 < 0i < 0.5. 
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1 INTRODUCTION 

Measurements of peculiar motions provide a fundamental tool to probe the mass distribution in the local universe. In 
the linear regime of gravitational instability, a simple relation between peculiar velocity, v, and mass density contrast, 
5m, can be easily obtained from mass conservation, either in differential 



where Q m is the cosmological mass density parameter. Together with the commonly used simplifying assumption 
of linear biasing, S g — b g 5 m , where 8 g is the galaxy density contrast and b g the galaxy biasing parameter, these 
two equations provide a relation between observable quantities: the peculiar velocity of luminous objects, v, and the 
density contrast, 5 g , which can be obtained from large all-sky redshift surveys. The IRAS Point Source Catalogue 
(Beichman et al. 1998) is particularly suited for this purpose because of the good sky coverage and homogeneity of 
the survey. In this paper, we present a comparison of peculiar velocity data with a redshift survey based on the IRAS 
catalogue. Throughout the paper, we will use the subscript j to indicate that an IRAS density field has been used. 

Comparing 5 g with v using any of the two equations above allows us to estimate the quantity [3 g = Q^' /b g and to 
test the validity of the gravitational instability hypothesis. Although the two equations (Q) and (||) are mathematically 
equivalent, they lead to two different strategies for measuring j3 g . Equation [j] is used to perform the so-called density- 
density (d-d) comparisons which typically consist of the following steps: a 3-D velocity field reconstruction from 
observed radial velocities; differentiation of v(r) and use of eqn. (Q) to compute 8 g ; comparison to the observed galaxy 
density fields. The first step is the least trivial and requires some additional theoretical assumptions. Bertschinger 
& Dekel (1989) successfully implemented a d-d technique by assuming that v(r) is irrotational, in the POTENT 
reconstruction method. The many applications of POTENT to various datasets have consistently led to large values 
of pi consistent with unity (see Sigad et al. 1998 and references therein). Equation (Q) is at the core of the so-called 
velocity velocity (v-v) comparisons. In this approach, one computes the mass density field obtained from the galaxy 
distribution in the redshift survey, uses equation ^ to predict a peculiar velocity field and then compares it with 
the observed galaxy velocities. The v-v methods have been applied to most of the catalogues presently available and 
have given values of fli which are typically in the range 0.4 — 0.6 (see Willick 2000, for an updated summary of the 
various results). 

The v-v methods are commonly regarded as more reliable and robust than the d-d ones because they require less 
manipulation of the data. The values of f3i obtained by these analyses are significantly smaller than unity, irrespective 
of the velocity tracers, model gravity field and comparison technique used. However, some of the v-v analyses showed 
evidence for a poor match between models and data, which would render the estimate of j3i meaningless. Davis, Nusser 
& Willick (1996) used their ITF technique (Nusser & Davis 1994) to compare the gravity field derived from the IRAS 
1.2 Jy survey (Fisher et al. 1995) with the peculiar velocities obtained from the Mark III catalogue (Willick et al. 
1995, 1996, 1997a). The coherent dipole residuals that they found were taken as evidence for significant discrepancies 
between modeled and observed velocity fields. Willick et al. (1997b, VI hereafter) and Willick & Strauss (1998, 
V2 hereafter) also considered the IRAS 1.2 Jy velocity predictions and the Mark III dataset but compared them 
using the VELMOD method. They were able to obtain a good fit to the data only by introducing a physically 
motivated, external quadrupole contribution to the model velocity field. Da Costa et al. 1998 (D98 hereafter) found 
good agreement between the peculiar velocities of galaxies in the SFI catalogue (Haynes et al. 1999a, 1999b) and 
those derived from the IRAS 1.2Jy gravity field, by performing an ITF comparison. The same ITF method has been 
recently applied to compare two different datasets: the IRAS PSCz redshift survey (Saunders et al. 2000) and the 
peculiar velocities in the recently completed ENEAR catalogue (da Costa et al. 2000). Also in this case the agreement 
between model and data was satisfactory (Nusser et al. 2000). 



V ■ v = -Q° m b S, 



(1) 



or integral form, 




(2) 
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In this paper, we use the VELMOD technique to compare the PSCz velocity prediction to the SFI dataset. As 
for any v-v comparison our main aim is to constrain j3i and to investigate the adequacy of the PSCz model velocity 
field. However, rather than simply adding one more measurement of j3i to those already in the literature from other 
v-v comparisons, we hope to address some more specific questions which should help simplify the rather complicated 
picture that has emerged from the results of the various v-v comparisons. Our goal is to check whether the PSCz- 
IRAS gravity field still provides a good fit to SFI data when the comparison is performed with the VELMOD method 
rather than the ITF. We also want to exploit fully the denser and deeper PSCz catalogue and see whether we can 
improve the agreement between the gravity field and the measured velocities, thus reducing the uncertainties in the 
estimate of /3j. 

In Section || we review the basics of the VELMOD technique and describe its current implementation. The 
SFI sample is presented in Section ^| and the IRAS PSCz catalogue and its gravity field are described in Section ^. 
VELMOD is tested in Section ^ and the results of its application to the SFI catalogue are presented in Section |(S| An 
analysis of errors based on the magnitude and velocity residuals is performed in Section Finally, in Sections ^ and 
^| we discuss the results and conclude. 



2 THE VELMOD METHOD 

VELMOD is a maximum likelihood method introduced by VI and described in detail in VI and V2. Here we 
simply outline the main points of the method, focusing on its implementation in the case of a forward TF relation. 
In Strauss & Willick's (1995) terminology, VELMOD uses a Method II approach, i.e. takes the TF observables 
(apparent magnitude and velocity width) and the redshift of an object and quantifies the probability of observing 
the former given the latter, for a particular model of the velocity field and a TF relation. This probability is then 
maximized with respect to the free parameters of the velocity model and the TF relation. Unlike previous Method 
II implementations (e.g. Hudson 1994), VELMOD replaces the redshift-distance relation with a joint probability 
distribution of distance and redshift. This probabilistic approach allows a statistical treatment of all those effects 
(small-scale velocity noise, inaccuracy of the velocity model and existence of triple-valued regions) that spoil the 
uniqueness of the redshift-distance mapping. 

VELMOD does not require smoothing of the TF data which, along with the allowance for triple-valued regions 
and small-scale velocity noise, allows one to probe the velocity field in high-density regions, thus exploiting the denser 
sampling of the new PSCz galaxy catalogue. Another convenient feature of VELMOD is that it does not require an 
a priori calibration of the TF relation, which is a common issue of concern in peculiar velocity studies. Instead, a fit 
of the parameters of the TF relation is performed simultaneously with the parameters of the velocity field. 

2.1 Implementation of VELMOD 

For each galaxy, the angular position (l,b), redshift (cz), apparent magnitude (m), and velocity width parameter 
(rj = log 10 (W) — 2.5, where W is twice the rotation velocity of the galaxy) are taken from the SFI catalogue. 
Hereafter, we will always use the redshift measured in the Local Group (hereafter LG) frame unless otherwise specified. 
In addition, a model for the density field is needed to account for inhomogeneous Malmquist bias. Such a model is 
obtained from the distribution of IRAS PSCz galaxies, together with a self-consistent model for the peculiar velocity 
field (as described in Section ^). The velocity model is completely specified by the value of the f3i parameter and 
a one- dimensional velocity dispersion, o~ v , which quantifies both the inaccuracy of the velocity model and the true 
velocity noise arising from small-scale nonlinear motions. Strauss, Ostriker & Cen (1998) and V2 have shown that the 
velocity dispersion on small scales is an increasing function of the local density. The SFI galaxy sample considered 
in this work consists of field spirals which avoid high-density environments and thus we ignore such dependency and 
assume a constant a v , independent of the environment. 

The probability of observing a magnitude m for a galaxy with a given r] and distance r is modeled by a linear 
TF relation, 
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m = M{rj) +51og(r) = A TF - b TF r] + 51og(r) + a r , (3) 

where o> is a Gaussian random distribution with zero mean and dispersion otf- This TF relation is completely 
specified by its zero point (Atf), slope (&tf) and scatter (cttf). In this work we do not model a possible dependency 
of the TF scatter on the luminosity or velocity width. In section 6.4, we will find that our results are insensitive to 
this approximation. 

The final ingredient needed to compute the probability of observing the measurements of the SFI catalogue are 
the selection function of the observational quantities and their correlations. Systematic errors can affect peculiar 
velocities computed with the forward TF relation if selection effects are not propery accounted for. For that purpose, 
we have used the correlation and selection models of Freudling et al. (1995). 

With the above assumptions, one can compute the probability that the i-th object of the sample with recession 
velocity czlg and velocity width parameter 77 will have an apparent magnitude m: Pi(m\r), czlg). To evaluate this 
conditional probability one needs to integrate the joint probability distribution Pi(m,r), czlg) over m. Although 
analytic approximations for this integral have been introduced by V2, which are valid away from triple-valued regions 
and at distances much larger than a v , in this work we perform an explicit numerical integration for all galaxies. 

From the overall probability, which we obtain by multiplying the single-object probabilities P = J~| Pi(m\r), czlg), 

we compute the likelihood C = — 21nP which is then minimized at each of 19 values of /3j = 0.1, 0.15 0.95, 1.0 by 

continuously varying the remaining free parameters (i.e. the three TF parameters and a v ). The function £. m in(/3i) 
obtained by this procedure is then fitted with a cubic function and the maximum- likelihood value of /3r, f3 m in, is 
found at the minimum of the curve. Extensive tests with mock catalogues performed by VI have shown that /3 m m is 
an unbiased estimator of the true fii parameter. 

VI and V2 found systematic residuals between the IRAS 1.2Jy predicted and the Mark III observed velocities 
which they modeled as a velocity quadrupole with a distance-dependent amplitude. This quadrupole, which they 
treated as a free parameter in their VELMOD analysis, is likely to arise from the missing contributions to the model 
velocity field from the mass distribution in the regions beyond the limits of the IRAS 1.2Jy survey and from shot 
noise. The PSCz survey should be deep and dense enough to reduce such discrepancies allowing us to drop this 
extra parameter. Therefore, we choose not to allow for any external contribution to the PSCz velocity field in our 
VELMOD analysis. In Section ^ we will show that the PSCz velocity model within 6000 kms -1 does indeed constitute 
an acceptable fit to the SFI peculiar velocity field, with no need to introduce external contributions. Finally, unlike 
VI and V2, we do not model the uncertainties in the LG velocity as a free parameter. Instead, we shall test explicitly 
the robustness of our results by running a few VELMOD experiments in which the LG velocity is allowed to vary 
within its 1-a error range. 



3 THE SFI SAMPLE 

The SFI sample of galaxies (Giovanelli et al. 1997a; Haynes et al. 1999a,b) is a homogeneous all-sky sample of 
galaxies for which I-band Tully-Fisher parameters are available. The sample uses new observations for declinations 
S > —40° and re-reduced data from Mathewson, Ford and Buchorn (1992) in the southern polar cap. The sample is, 
by design, angular diameter limited, with different limiting diameters, Du m , for different redshift shells. The limiting 
diameters, expressed in unit of 0.1 arcmin, are Du m = 25, 16 and 13 in the redshift ranges czlg < 3000 kms -1 , 
3000 < czlg < 5000 kms -1 and 5000 < czlg < 7500 kms -1 , respectively. However, for a number of different reasons, 
some bright galaxies smaller than the stated diameter limit were included in the sample. These additional galaxies, 
which amount to ~ 15% of the total, do not constitute a strictly magnitude-limited sample. All of them are brighter 
than an apparent Zwicky magnitude of m 2 = 14.5. In a detailed investigation of these extra galaxies, we have not 
found any dependence of the selection on distance or apparent diameter. This is different from the distance-dependent 
selection function that Willick et al. (1996) have used to describe the MAT sample. For the present analysis, we have 
therefore approximated the SFI selection sample as the combination of a strictly diameter-limited sample with the 
above criteria, and a magnitude-limited sample with limiting Zwicky magnitude of 14.5. For such a case, Willick 
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(1994) has derived the VELMOD formalism. In his terminology, this is the Two-Catalogue Selection case and the 
selection function, expressed in a form suitable for the VELMOD analysis, is given by eq. (63) of Willick (1994). We 
have found that the results presented in this paper do not change significantly when using m z — 13.5 for the putative 
magnitude limit of the extra bright galaxies 

Freudling et al. (1995) have modeled the correlations between the quantities used to define the selection criteria, 
(m z ,D), and the TF observables, (mi,rj), as 

m z — 1.6 + mi + 0.5-q with dispersion <r m _ = 0.59 (4) 
and 

D = 3.82 + -0.205m/ - 0.1027? with dispersion a D = 0.121, (5) 

where mi indicates the I-band apparent magnitude. These correlations can be easily translated into the VELMOD 
formalism. 

Below, we will consider the following sub-samples drawn from the SFI catalogue. Most of the analysis is carried 
out with the 989 galaxies with czlg < 6000 kms -1 and r\ > —0.25. These cuts allow us to perform a homogeneous 
comparison with the ITF analysis of D98. In addition, we will consider separate subsamples restricted to three 
independent redshift-space shells 2000 kms -1 thick and with external radii of 2000, 4000 and 6000 kms -1 . These 
subsamples mix different nominal selection criteria and serve as an additional check that selection effects are well 
accounted for. The three subsamples contain 158, 355 and 496 galaxies, respectively. 

SFI galaxies avoid rich clusters and high density environments. In particular none of them belong to the Virgo 
cluster. Therefore, there is no need to adopt any grouping procedure and the hypothesis of a small-scale velocity 
dispersion, o~ v , independent of the environment is well justified. 



4 MODEL DENSITY AND VELOCITY FIELDS FROM THE PSCZ SURVEY 

The models for the density and velocity fields are obtained from the distribution of IRAS galaxies in the recently 
completed PSOz all-sky redshift survey (Saunders et al. 2000). This catalogue, which basically extends the old 1.2 Jy 
one (Fisher et al. 1995), contains ~ 15500 IRAS PSC galaxies with a flux at 60 /im larger than 0.6 Jy. A complete 
description of the dataset, its selection criteria, the procedures adopted to avoid stellar contamination and galactic 
cirrus, are given in Saunders et al. (2000). For our purposes, the most interesting features of the catalogue are the 
large area sampled (~ 84% of the sky), its depth (the median redshift is 8500 kms -1 ) and dense sampling (the mean 
galaxy separation at 10000 kms -1 is (I) ~ 1000 kms -1 j^j compared with a value of (l) ~ 1500 kms -1 for the IRAS 
1.2 Jy catalogue). 

The flux-limited nature of the catalogue causes the number of objects to decrease with distance. This decrease 
is quantified by a radial selection function. Various authors have used different estimators to compute the selection 
function (Springel 1996, Canavezes et al. 1998, Sutherland et al. 1999, Branchini et al. 1999, Monaco and Efstathiou 
1999, Saunders et al. 2000). Branchini et al. (1999, hereafter B99) found that different selection functions induce 
variations smaller than 5% in the model density and velocity fields within a distance of 20000 kms -1 . In this work 
we use the selection function specified by eqn. (1) of B99. 

The density and velocity field in real space are obtained from the redshift space distribution of PSCz galaxies by 
implementing Method 1 of B99, which uses the iterative technique of Yahil et al. (1991) to minimize redshift-space 
distortions. The procedure relies on gravitational instability theory, assumes linear biasing, and is valid in the limit of 
small density fluctuations where linear theory applies. At each step of the iteration the gravity field, g, is computed 
from the distribution of the 11206 PSCz galaxies within 20000 kms -1 . The gravity field is subsequently smoothed 
with a top hat filter of radius ~ 500 kms -1 which allows us to assume linear theory and thus to obtain the smoothed 

* Throughout this paper we measure distances in velocity units (kms -1 ). This is equivalent to setting the Hubble constant 
equal to unity. 
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peculiar velocity of each galaxy from the acceleration, v oc /3/g. The new distances of the objects, r, are then assigned 
assuming a unique redshift -distance relation r — cz — u, where u is the radial component of the peculiar velocity 
vector. The procedure is repeated until convergence is reached. 

The final products are the real space positions of the galaxies and their velocities for a given value of /3i. We 

have ran 19 different reconstructions with /3i = 0.1, 0.15 1.0. The continuous density field is obtained by smoothing 

the galaxy distribution on 129 3 points of a cubic grid inside a box of 19200 kms" 1 a side with the LG at the centre, 
using a Gaussian filter of 300 kms -1 (G3, hereafter). The associated smoothed velocity fields are computed for the 
appropriate value of /3j from eqn. ([!]). 

The procedure returns 19 models of the density field (one for each value of f3i), along with their associated 
peculiar velocity fields, both defined at the gridpoint positions. Velocity predictions are made in the LG frame to 
minimize the uncertainties derived from the lack of information about the mass distributions on scales larger than 
the size of the PSCz sample. 

5 TESTING THE VELMOD IMPLEMENTATION 

In this work we rely on the error analysis performed by VI, based on an extensive application of VELMOD to realistic 
mock catalogues, which we do not repeat here. Instead, we present the results of two tests aimed at checking the 
reliability of our VELMOD implementation. 

5.1 Tests with Ideal PSCz and SFI Mock Catalogues 

As a first test, we applied VELMOD to a suite of SFI and PSCz mock catalogues. Errors in the VELMOD analysis are 
dominated by inhomogeneous Malmquist bias, with uncertainties in the model velocity fields playing an important 
role only in the innermost part of the sample. To account properly for Malmquist bias, the mass distribution in the 
mock catalogues needs to mimic the one in our local universe as close as possible. Unlike VI, however, we do not 
extract mock catalogues from constrained N-body simulations (i.e. similar to those produced by Kolatt el al. 1996). 
Instead, we construct 'ideal' PSCz and SFI mock catalogues from the G3 smoothed PSCz density and reconstructed 
linear velocity fields (described in the previous section) for fii — 1.0. Nonlinear effects were mimicked by adding a 
constant Gaussian noise of 150 kms -1 to each of the Cartesian components of the model velocity vectors. Such a 
small value is meant to mimic the observed "coldness" of the cosmic velocity field (e.g. Strauss, Ostriker and Cen 
1998). 

Twenty mock PSCz catalogues where generated with Montecarlo techniques using different random seeds. Se- 
lection criteria close to the observational ones were applied to mimic the PSCz selection function and the presence 
of unobserved regions. Mock PSCz galaxies were generated assuming a probability proportional to the local density, 
i.e. assuming that they trace the mass distribution (bi = 1). The same reconstruction method used for the real 
data was then applied to each of these mocks to obtain 20 mock PSCz velocity models for each of the 20 values of 
Pi = 0.5,0.55 1.5. 

Similarly, we generate 20 mock SFI catalogues by Montecarlo resampling the original PSCz G3-smoothed density 
and velocity fields with the added thermal noise. First, catalogues with large numbers of galaxies were created. Next, 
absolute magnitudes were assigned according to the luminosity function. Subsequently, diameters and velocity widths 
were assigned according to the correlations described in Section ^| TF parameters close to the ones found in our final 
analysis were used to assign velocity widths. Finally, the nominal selection criteria of the SFI sample were used to 
select galaxies for the final catalogues. The additional galaxies mentioned in Section |§| were simulated by randomly 
including galaxies up to the magnitude limit until the number of galaxies not satisfying the nominal selection criteria 
matched the one of the SFI sample. The resulting mock samples nicely reproduce the observed redshift distribution 
of the SFI sample. 

The results of applying VELMOD to the 20 SFI mock catalogues are summarized in table [I]. The true values of 
the free parameters used in the mock catalogues are listed in the first row while those obtained from the VELMOD 
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analysis are shown in the second row. For each parameter we report the average value from the 20 mocks, the error in 
the mean and, in parenthesis, the typical error in a single realization. These results indicate that VELMOD returns 
an unbiased estimate of the free parameters. 

Given the ideal nature of the velocity field in these mock catalogues the errors displayed in table |l| are likely 
to underestimate the real ones. In what follows, however, we will be mainly interested in the errors on (3i which, as 
VI demonstrated, can be obtained from the values at which Cmin(Pi) differs by one unit from its minimum value at 

f3min • 

5.2 Tests with the Mark III Catalogue 

In the second test, we repeated part of the VELMOD analysis performed by VI and V2, but using the PSCz model 
velocity field instead of the IRAS 1.2 Jy one. We applied VELMOD to the Aaronson et al. (1982, hereafter A82) 
and Mathewson et al. (1992, MAT) subsamples selected according to the VI prescriptions. The selection functions 
for A82 and MAT and their coefficients were taken from Willick et al. (1996). Both samples are spatially limited to 
a redshift czlg = 3000 kms -1 . A third subsample we use, MAT2, also obtained from the Mathewson et al. (1992) 
catalogue, coincides with the one considered by V2 and extends out to czlg = 7500 kms -1 . 

The results are summarized in table ^ where the values of the free parameters obtained from our VELMOD 
analysis are compared to those obtained by VI and V2 (in parenthesis). Only errors in fii are quoted. Our results 
generally agree with those of VI and V2. The analysis of the MAT and MAT2 samples returns values of /3i smaller 
than those obtained by VI and V2. The difference is within the la error bar and probably reflects the difference 
between the IRAS 1.2 Jy and PSCz model velocity fields. 

The largest discrepancy is the a v parameter for which we consistently obtain values larger by 30 — 70% than those 
estimated by VI and V2. Such discrepancy is statistically significant since the expected errors on a v estimated from 
the mock catalogues and, more qualitatively, from the likelihood curve, are of the order of 40 kms . To interpret 
this discrepancy one has to keep in mind that the value of a v is approximately given by the sum in quadrature of 
random velocity noise (or) and the uncertainties in the model velocity field (o\m). 

The only possible explanation for the discrepancy in a v is that or, which is an intrinsic property of the velocity 
catalog, is the same in both analyses, but the PSCz velocity model has a larger gm- Indeed, since B99 estimate 
<tm ~ 130 kms -1 for the PSCz model, while VI find <tm ~ 84 kms -1 for the IRAS 1.2Jy model, one may think that 
the intrinsic velocity noise would then be similar in the PSCz (or ~ 120 kms -1 ) and IRAS 1.2Jy (ot ~ 100 kms -1 ). 
However, it is difficult to understand why errors in the PSCz model should be larger than in the IRAS 1.2Jy one. 
The denser sampling and larger depth of the PSCz catalogue should decrease the uncertainties in the model velocity 
field, not increase them. In fact, the errors quoted by VI and B99 were estimated using different mock catalogues. 
VI used mocks derived from the reconstruction technique of Kolatt et al. (1996) which are based on a PM N-body 
code and are intrinsically "colder " than those used by B99 which were taken from the higher resolution AP 3 M 
N-body simulations of Cole et al. (1998). Since the precision of the velocity reconstruction decreases with increasing 
small-scale random velocity noise, it is reasonable to suspect that the difference in the values of cta/ estimated by VI 
and B99 simply reflects the differences in the mock catalogues rather than differences in the errors of the two model 
velocity fields, and that an error analysis based on the same set of mock catalogs would show that the value of ctm 
for the PSCz model is comparable, if not smaller, than for the IRAS 1.2Jy model velocity field. 

Nevertheless, the ct„ discrepancy must reflect some difference in the IRAS PSCz and 1.2Jy velocity prediction. 
The only remaining possibility is that such a difference is systematic rather than purely random. Indeed, the two 
models are different because they are based on two different redshift catalogues and because they are derived using 
two different reconstruction methods. The PSCz catalogue allows a denser sampling and therefore traces small-scale 
density fluctuations which are washed out by shot noise in the IRAS 1.2Jy catalogue. The IRAS 1.2 Jy model used 
by VI and V2 was obtained using the technique of Sigad et al. (1998) which differs from the one of B99 in the 
treatment of triple-valued regions, the use of Wiener filtering, and, most importantly, in the allowance for mildly 
nonlinear motions. As a result, the PSCz velocity model is intrinsically more linear than the IRAS 1.2Jy one which 
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Figure 1. The VELMOD likelihood , C m i„(0i), (filled symbols) for the full SFI sample. The dotted curve represents a cubic 
fit to the points. 

also lacks power on small scales. The linearity of the PSCz model forces any discrepancy between true and model 
velocities to contribute to a v . In the IRAS 1.2Jy case, such discrepancies are, on average, smaller because of the lack 
of small-scale power and because they are partially absorbed by the mildly nonlinear motions. The net result is that 
VELMOD analyses based on the IRAS 1.2 Jy velocity model return values of a v systematically smaller than those 
based on the PSCz velocity model and a value of /3 slightly larger. 



6 RESULTS 

In this section we present the results of applying the VELMOD analysis to the SFI sample within 6000 kms -1 and 
we check the robustness of our results using smaller subsamples, smoothing scale and errors in the LG velocity. 



6.1 Application to the Full SFI Catalogue 

We ran VELMOD on the full SFI sample of 989 galaxies with redshift czlg < 6000 kms^ 1 and 77 > -0.25. The 
PSCz density and velocity fields used in this calculation were smoothed with a Gaussian filter of effective radius 300 
kms -1 . The likelihood was minimized at each of 19 values of f3i = 0.1,0.15, ....1.0 by varying the free parameters 
Atf, brF, &tf and a v . The resulting values of £min(/3i) are represented by filled dots in fig [jj In order to find the 
actual minimum, we fit a cubic function to C m in{f3i)- 

Cmin = C-0 + dl(/3j — Pmin) 2 + Cb^^Pl ~ Pmin) 3 (6) 

The value of (3 m in and its errors are indicated in the plot and are also listed in table |^ together with the values of 
the other parameters found at the minimum of the likelihood curve. 
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Figure 2. Best fit parameters obtained from running VELMOD on the full sample as a function of TF slope (top left), TF 
zero point (bottom left) velocity dispersion (top right) and TF scatter (bottom right). 



The values for the TF parameters listed in table ^ should be compared with the values found by Giovanelli et al. 
(1997b, G97) for an associated sample of cluster galaxies, termed SCI, selected in a similar way to the SFI sample, 
which are commonly used as the assumed TF relation for the SFI sample. They find zero point, Atf, of-6.00±0.02Q 
and slope, 6tf, of 7.5 ± 0.2 with a scatter of about 0.36 magnitudes. While the range of values for the slope of this 
determination from cluster galaxies overlaps with the current results, we obtain a smaller value for the zero point and 
a larger scatter. This is in agreement with other analyses of the SFI dataset (D98, Freudling et al. 1999). A similar 
difference in the TF scatter of cluster and field samples has been found in other datasets (e.g. Bothun & Mould 1987; 
Freudling, Martel & Haynes 1991). 

The recovered value of a v = 250 kms -1 is twice as large as the value obtained by VI and V2. As we have already 
discussed, part of this discrepancy may simply reflect differences between our model velocity field and the one used by 
VI and V2. It may also indicate that the PSCz velocity model does not provide a satisfactory fit to the SFI dataset. 
We will return to this important point in Section 0. 

Figure ^ shows the variation of the four parameters with j3i, where the other parameters are allowed to vary. 
The velocity dispersion reaches a minimum when /3i ~ f3min, a behavior consistent with the analysis of the mock 
catalogues performed by VI (but not with their analysis of the real A82+MAT sample). The TF parameters are very 
insensitive to /3j. The largest variation (~ 5%) occurs for the TF scatter which has a minimum around /3i = 0.55, 
significantly larger than j3 m i n . This is not alarming since, as stressed by VI, minimizing the TF scatter does not 
correspond to maximizing the likelihood. This is probably due to the covariance between a v and cttf: the increase 
of a v after its minimum is compensated by a further decrease of cjtf- 



t In this work absolute magnitudes refer to distances in units of kms 1 , and not of 10 pc, as often adopted in the literature 
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6.2 The Effect of Smoothing 

In their analysis, VI noticed that smoothing the model velocity field with a G5 filter had little effect on the VELMOD 
results. Along with their small value of a v , this was taken as an indication that density fluctuations on scales between 
300 and 500 kms -1 contribute little to both the Mark III and the IRAS 1.2Jy velocity fields. However, the larger value 
of Ov we have found in our VELMOD analysis of the A82, MAT and SFI samples seems to indicate that contributions 
from small scales are non-negligible. 

To investigate this issue, we have repeated the same exercise and performed a VELMOD analysis using a G5 
smoothed PSCz model velocity field instead of the G3 one. We find that the results are insensitive to the smoothing 
scale, with /3 min increasing only by ~ 5%, consistent at the la level with the value found in the G3 case. The same 
is true for the other free parameters with the exception of a v which increases by ~ 15%. Moreover, the value of 
Cmin at Pmin increases by 10 units with respect to the G3 case, corresponding to a probability decrease by a factor 
c 5 ~ 90 (as C — — 21nP). How do we interpret these results ? With a G5 filter, the model velocity field does not 
receive contributions from scales smaller than ~ 500 kms - and larger values of f3i and a v are needed to match 
the amplitude of the observed velocities. This means that the linear PSCz model velocity field does receive a non- 
negligible contribution from small scales. The significant increase in the likelihood indicates that the G3 velocity 
field provides a better fit to the SFI velocities, i.e. small-scale fluctuations do contribute to the peculiar velocities of 
galaxies in the real world. 

6.3 Uncertainties in the Velocity of the Local Group 

Since model predictions are given in the LG frame (i.e. the LG velocity is subtracted from all other velocities), we need 
to quantify the impact of uncertainties in the predicted and measured LG velocity when performing the VELMOD 
analysis. VI and V2 tackled the problem by introducing a free parameter to model a random component to the 
LG velocity vector. Here we take a different approach aimed at minimizing the number of free parameters in the 
analysis. Instead of modeling the uncertainties in the LG velocities, we quantify the impact that these uncertainties, 
independently estimated, may have on our VELMOD analysis. Errors in the LG velocity derive from uncertainties in 
modeling the LG velocity and in transforming redshifts from the heliocentric to the LG frame. The former have been 
quantified by Schmoldt et al. (1999), while for the latter we consider the recent work of Courteau & van den Bergh 
(1999). We assume that these error estimates are independent and so we compute the total uncertainty by adding 
them in quadrature. The resulting total error in the LG velocity is ~120/3j kms -1 for each Cartesian component. 

To evaluate the impact of these errors we ran 10 VELMOD analyses in which the LG velocity was perturbed 
with a Gaussian noise of the same amplitude. The effect is to increase the random errors without introducing any 
systematic bias. Averaging over the 10 experiments we obtain a value of /3j = 0.43 with a 1-a error around the 
mean of 0.05, slightly larger than the typical error in a single realization. Similar considerations apply to the other 
parameters. In all but one experiment the value of jCmin(f3min) was larger than in the unperturbed case. In the one 
case with smaller likelihood (by only two units) the perturbation turned out to be very small, in agreement with the 
VI and V2 results in which the extra random components added to the LG velocity vector turned out to be trivially 
small. 

6.4 Breakdown by Rotation Velocities 

In their calibration of the TF relation for SCI galaxies, G97 found that the TF scatter is a function of the velocity 
width: rapidly rotating galaxies have a smaller TF scatter than slow rotators. A similar trend has also been detected 
by Federspiel, Sandage & Tamman (1994) and by Willick et al. (1997a) in some of the Mark III subsamples. In the 
SFI catalogue itself this effect has been included for the direct TF relation by Freudling et al. (1999), but not in 
the inverse TF relation (D98). In the current analysis, we have chosen to neglect this effect. To estimate the impact 
of this approximation we divided the SFI catalogue into two subsamples according to the rotational velocity of the 
galaxies. In the first catalogue, we included only objects with a velocity width parameter r\ > —0.1, while in the 
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Figure 3. Likelihood curves, £ m i n (/3i), for three different redshift intervals. In each case the value of /3„ 
galaxies in the samples are indicated. The dotted curves are cubic fits to the C m i n (f3i) points. 



and the number of 



second we included slow rotators with r\ < —0.1. The results of applying VELMOD on these two samples are shown 
in table []. The TF dispersion is found to be smaller for fast rotators, as could be expected. However, the values of the 
other free parameters, in particular /3 m in, are very similar in the two subsamples, showing that our approximation of 
a constant otf has little impact on our /3i estimates. 



6.5 Breakdown by Redshift 

As a last robustness test we have divided our sample into three independent redshift shells 2000 kms -1 thick and 
applied VELMOD to each of them. Since the selection criteria are different in the three shells, this test serves as a 
check of whether selection effects are properly taken into account by our procedure. The resulting £ m i n (/3j) curves are 
displayed in fig. |^ and the free parameters are listed in table [| The values of (3 m in in the three shells are consistent 
with each other. The remaining parameters are also in good agreement, with only two exceptions. One is the TF slope 
which in the innermost shell is ~ 10% shallower than in the rest of the sample. The analysis of SFI mock catalogues 
reveals that such a discrepancy is significant at the 1.5c level. However, as we have already pointed 



in Section 5.1 



out, the errors estimated using those catalogues are smaller than the real ones. A VELMOD analysis of more realistic 
SFI mocks would return larger errors and decrease the statistical significance of the slope variation. 

The second, more serious, discrepancy is observed for er„ which in the outermost shell increases by a factor of 
~ 2.7. Such a variation is much larger than what we observe in the mock catalogues and deserves some further 
investigation. We have therefore sliced the SFI sample into thinner redshift shells, all of them 1000 kms -1 thick 
except the innermost one which we extended to 2000 kms -1 so that a comparable number of objects is contained 
in each shell. We then repeated the VELMOD analysis in each of these shells. The results are visualized in the 
histograms of fig ^. The two dashed lines represent the la error about the mean. The upper plot shows the radial 
behaviour of a v . Within 4000 kms -1 , a v ~ 200 kms -1 , comparable to what we have obtained from the analysis of the 
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Figure 4. Dependency of u v (upper plot) and m in (lower plot) on redshift. The histograms show the results of running 
VELMOD on 6 different redshift shells. In both plots the dashed lines represent lc errors around the mean. 

A82, MAT and MAT2 samples. This value is larger than those found by VI and V2 and, as discussed in Section |5.2] , 
this is probably because the the PSCz model velocity field is more sensitive to power on small scales. Around 5000 
kms -1 <j v increases to ~ 600 kms -1 and then stabilizes at a value of ~ 400 kms -1 . 

In order to test whether any particular part of the sky is responsible for this dependence of a v on distance, we 
have cut the sample into several complementary hemispheres (e.g. above and below the the Galactic plane, above and 
below the Supergalactic plane and so on) and, for each pair, we have ran VELMOD on the two sets of redshift shells. 
In all cases we have found that a v increases beyond 4000 kms -1 , suggesting that no particular cosmic structure is 
responsible for the increase in a v . 

The ability of VELMOD to constrain a v is a rapidly decreasing function of the redshift and we need to assess 
the statistical significance of the jump at 4000 kms -1 , i.e. we need to estimate the errors on a v . We compute them 
from the distribution of a v found from the analysis of our 20 mock SFI samples. This test reveals that the increase 
of c v is significant at about the 3cr level, leaving little doubt about the reality of the change at around 4000 kms -1 . 
However, as we have already stressed, such tests are only indicative since non-linearities are not properly modeled in 
our mock samples. 

To further assess the reality of this peak we performed a second test. We ran two VELMOD analyses using 
only objects with redshifts in the range 4000-6000 kms -1 . In the two calculations all the parameters were fixed to 
their maximum likelihood values except a v which was set equal to 200 kms -1 in one case and to 535 kms -1 (i.e. 
its maximum likelihood value) in the other. The smaller a v results in an increase of 12 units in the likelihood or a 
decrease in probability of ~ e 6 ~ 400, indicating that the increase in cr v seen in the external shell is indeed significant. 
It is therefore possible that the dramatic increase in a v reflects a disagreement between the velocity model and the 
SFI data on large scales. We will investigate this possibility further in the next section. 

The lower plot of fig ^ shows how p m in varies with redshift. VELMOD returns a very robust estimate of (3 m i„ 
which does not change significantly even when a v increases. Despite this encouraging evidence, it is sensible to 
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Figure 5. The VELMOD likelihood, C min (/3i), (filled symbols) for SFI galaxies with cz LG < 4000 kms" 1 . The dotted curve 
represents a cubic fit to the points. 



adopt a more conservative approach and repeat the VELMOD analysis within 4000 kms -1 , i.e. in the volume where 
there are no obvious indications of a possible mismatch between model and data. The results, listed in table [j, 
are in good agreement with those of the full sample. In particular, the minimum of the likelihood curve (shown in 
fig. |^), p m i n — 0.45 ± 0.05, is fully consistent with the results from the other subsamples. The velocity dispersion is 
~ 200 kms" 1 , in good agreement with those obtained from the analysis of the A82, MAT and MAT2 samples. 



7 ANALYSIS OF THE VELOCITY AND MAGNITUDE RESIDUALS 

VELMOD is a maximum likelihood technique in which the sources of variance, a v and (Ttf, are treated as free 
parameters. For this reason, the VELMOD analysis can only tell us which are the best values of j3i, Atf, 6tf, 
<ttf and a v for a given velocity field model, but it cannot address the question of whether the velocity model is an 
acceptable fit to the data. In this respect, the increase of a v at large radii found in the preceeding section is only 
suggestive of a mismatch between observed and modeled velocities, and a proper error analysis is needed to assess its 
statistical significance. In this section we address this problem by inspecting the magnitude and velocity residuals, 
following the formalism and notation of VI. 

7.1 Maps of Velocity Residuals 

We define the normalized magnitude residuals for each object of magnitude m: 
_ m ~ E(m\ri,cz) 

where the expected apparent magnitude, E(m\rj, cz), and the dispersion around it, A m , can be obtained by integrating 
over the constrained probability function, P(m\r), czlg), which is computed in the VELMOD analysis. The magnitude 
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Figure 6. The sky projection in Galactic coordinates, as seen in the LG frame, of the smoothed VELMOD velocity residuals, 
U SFI — upsCz, for /3r = 0.4. Open circles indicate objects that are inflowing relative to PSCz velocity predictions; crosses denote 
objects that are outflowing. The size of the symbols is proportional to the amplitude of the velocity vector. 

residuals where smoothed on a scale S = d/5, where d is the PSCz predicted distance of the generic object, and then 
converted into smoothed peculiar velocity residuals, usfi — u-PSCz, according to the VI prescriptions. The resulting 
maps of the velocity residuals are shown in figs. [], ^ and |^ for three different values of /3i and in three different 
redshift shells. When looking at the maps one should bear in mind that coherence up to scales of ~ 35° is induced by 
the adopted smoothing and therefore systematic mismatches between the model and reality can only be revealed by 
coherence on much larger angular scales. The visual inspection of the maps clearly shows that the velocity residuals 
for j3i = 0.4 are less coherent and have smaller amplitudes than those of the models with 0i — 0.1 and /3i = 1.0. The 
residuals in the f3i = 0.1 map show a pattern similar to (3i = 0.4 case, but have a larger amplitude and coherence. 
When (3i — 1.0, the residuals grow even larger and the map exhibits a clear dipolar structure. The velocity residuals 
of the model with 0i = 0.4 look qualitatively similar to the ones between the SFI and IRAS 1.2Jy velocities for 
0i — 0.6, as seen in fig. 7 of D98. However, a one-to-one comparison between the two sets of maps is not possible 
since different smoothing procedures were used and since our residuals are computed at the PSCz predicted distances 
whereas the ones of D98 are computed at the redshift space positions. 

7.2 Residual Correlation Function 

A more quantitative assessment of the goodness of fit can be obtained by computing the correlation function of the 
unsmoothed magnitude residuals: 



1 




(8) 



N(t) 
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Figure 7. The sky projection of the smoothed VELMOD velocity residuals usfi — upsCz f° r Pi = 01. 




Figure 8. The sky projection of the smoothed VELMOD velocity residuals, ugpj — tipsczi f° r Pi = 1-0 
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Figure 9. The correlation function of magnitude residuals plotted for f3j = 0.1 (lower plot), /3j = 0.4 (middle plot) and f3i = 1.0 
(upper plot). The results refer to the full SFI sample. 

where JV(r) is the number of galaxy pairs with predicted separation ckj < T ± 100 km s -1 . This correlation function 
applies to normalized magnitude residuals, i.e., does not depend on a v and otf and is only sensitive to genuine 
correlation among residuals. 

In fig. ^ we show the correlation function of all SFI galaxies with czlg < 6000 kms -1 for the same three values of 
Pi used to produce the velocity residuals shown in figs. |^, Q and ^[ The error bars represent Poisson errors N(t)~ ' 5 . 
In the Pi — 0.4 model, the correlation function appears to be consistent with zero almost everywhere apart from 
some positive correlation at separations smaller than 500 kms , i.e. of the order of the smoothing scale of the model 
velocity field. A significant excess correlation on small and large scales is detected for Pi = 0.1 and, to an even greater 
degree, for Pi = TO. 

To quantify the goodness of the fits we compute the quantity 

N bins 

(9) 



2 \ " 

X« = 



^ N( n )' 



where £(t) = N(t)i(j(t) and Nbi ns is the number of independent distance bins in which ip(r) is computed. VI have 
shown that if the residuals are indeed uncorrelated on a scale r then £(t) is a Gaussian random variable with zero 
mean and variance N(t). Under this approximation, the x| statistic is distributed as \ 2 with number of degrees of 
freedom equal to the number of independent distance bins. Any correlation among residuals will result in a larger Xp . 
Extensive tests with mock catalogues performed by VI revealed that Xe indeed has properties similar to a \ 2 statistic, 
with the same variance, but with a mean of ~ 0.87 per degree of freedom rather than unity. We have computed the 
quantity Xs, f° r 10 values of Pi ranging from 0.1 to TO. The results are shown in fig [l^. The continuous, heavy line 
shows the expectation value for a \ 2 statistic with Nu ns = 60 degrees of freedom, while the dashed and long-dashed 
lines show, respectively, the la and 3a deviations from that value. The dotted line represents the expectation value 
of x| according to the VI correction, i.e assuming a number of degrees of freedom of 0.87 x 60 = 52.2. 
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Figure 10. The statistic x| plotted for various values of f3j. The heavy line shows the expectation value for a \ 2 statistic with 
the same number of degrees of freedom. The dashed and long-dashed lines show la and 3cr deviations from the expectation 
value. The dotted line shows the expectation value of x| corrected according to VI. 

Despite its limited discriminatory power, the x| statistic clearly indicates that PSOz velocity models with 0.3 < 
Pi < 0.5, and therefore also our best model according to the VELMOD analysis, /3i — 0.42 ± 0.04, provide an 
acceptable fit to the SFI velocities. All other models with smaller or larger values of f3i can be ruled out at a level 
> 3a. 

We are now in a position to address the question of whether the PSCz model velocity field provides an acceptable 
fit throughout the whole SFI sample, particularly in the external regions where a v is large. We do that by computing 
the residual correlation function in the same three redshift shells previously used. The results are displayed in fig. [ll] 
which shows the correlation functions, ip(r), for the three subsamples. 

The value of appears to be consistent with the null hypothesis of no correlation in nearly all the distance bins 
in each of the three redshift shells, with the exception of an excess correlation at small separation in the outermost 
shell. The value of the statistic x| m each shell is indicated in the plots. It is always smaller than the number of 
bins used in each shells and comparable to the expectation value for x| computed using the VI correction. These 
results show that, despite the large value of a v obtained from the VELMOD analysis, the residuals of the fit are not 
correlated, i.e. the differences in magnitude between model and data are randomly distributed. 



8 DISCUSSION 

The results of our VELMOD analysis generally agree with those of independent analyses, with the exception of the 
velocity noise, a v . We find a significantly larger value than that found by VI and V2 who compared the Mark III 
velocities to the IRAS 1.2Jy model velocity field. A VELMOD analysis of SFI subsamples limited in redshift revealed 
the existence of two regions. An inner sphere of radius 4000 kms -1 in which a v ~ 200 kms -1 , a value which is 
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Figure 11. Residual correlation function in three redshift shells for the model /3j = 0.4. The redshift ranges, values of xij 
number of bins and, in parenthesis, the number of degrees of freedom corrected according to VI are also shown in the plots 



believed to reflect the cumulative effect of velocity noise and errors in the model predictions (as corroborated by the 
results of the analysis of the A82 and MAT samples), and an external region in which a v ~ 500 kms -1 . 

If these large values of a v were caused by systematic deviations of the SFI velocity field from the model predictions, 
then the residuals should correlate on scales larger than the smoothing length. However, this is not the case (see 
fig. 10). In fact, the residuals turn out to be quite small and do not show any significant spatial correlation. This 
is true for all redshift shells, including the outermost one where cr v is large. Moreover, a quantitative analysis based 
on a \ 2 statistic shows that the PSCz velocity model with 0.3 < fli < 0.5 does provide a satisfactory fit to the 
SFI data without the need to introduce any external fields as in the VI and V2 VELMOD analyses. This can be 
interpreted as implying that discrepancies between the model and the data do exist, but not in the form of peculiar 
motions coherent on scales larger than ~ 300 kms . Rather, they are randomly distributed and thus properly 
quantified by means of the total variance, to which both a v and otf contribute. It is worth noticing that the increase 
in er„ occurs at the distance at which large structures like the Perseus-Pisces supercluster and the Great Attractor 
appear in the sample. Non-linear effects may be strongest in these structures, leading to a mismatch on small scales 
between measured and predicted velocity fields which could contribute to the increase in a v . To investigate whether 
these structures are indeed responsible for the increase in a v we have ran two VELMOD analyses in which we have 
excluded either SFI galaxies in the Perseus-Pisces supercluster (i.e. with 120° < I < 170°, —35° < b < 5° and 
4000 kms -1 < cz LG Ze6000kms -1 or in the Hydra Centaurus region (i.e. with 260° < / < 330°, -10° < b < 40° and 
2000 kms -1 < czlg < 6000 kms -1 ). A third analysis has been performed after excluding galaxies in both structures. 
The results turned out to be very similar to those shown in figure 4 for the full sample case (the differences are of the 
order of 20 kms -1 ). The abrupt increase of a v at czlg > 4000 kms -1 is present in all the cases explored, showing 
that neither the Perseus-Pisces supercluster nor the Great Attractor contribute appreciably to the variation of a v . 

A second possibility is that the increase in a v with redshift is caused by a comparable increase in the errors of 
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the velocity model. Indeed, B99 quantified such effect (see their eqn. 18) which, however, is far too small to explain 
the observed increase in a v . 

The final possibility is that the TF scatter decreases significantly with distance, which seems reasonable, since 
the scatter tends to decrease for high-linewidth galaxies, and these are seen preferentially at large distances. Because 
of the covariance between a v and otf (see VI), this could lead to an overestimation of a v and an underestimation 
of otf- In this context, it is interesting to note that a power spectrum analysis of the same SFI dataset, like the 
one performed by Freudling et al. (1999), also returns different results if computed separately for the closer and 
more distant halves of the sample. However, such analysis reveals that it is the inner part of the sample which has 
an apparently enhanced a v . Given these conflicting results, we consider it unlikely that the apparent increase in a v 
reflects a true property of the velocity field. In any case, it is worth stressing that the increase of a v with distance 
does not affect the estimate of /3/ which, as shown in fig. 4 and table 3, does not show significant variations. 

9 SUMMARY AND CONCLUSIONS 

We have implemented, tested and used the VELMOD technique (VI, V2) to compare the model velocity field obtained 
from the spatial distribution of PSOz galaxies with the velocities of field spiral galaxies within 6000 kms -1 in the SFI 
catalogue. This comparison allowed us to estimate the value of the Pi parameter, the amplitude of the small-scale 
velocity noise, a v , and to calibrate the TF relation of SFI galaxies. 

VELMOD returns an estimate of (3i which is very robust to various systematic and random errors that may 
enter the analysis at various stages. For the full sample, we obtain (5i = 0.42 ± 0.04, while a value of /?/ = 0.45 ± 0.05 
is obtained when a more conservative cut at 4000 kms -1 is applied to the SFI sample. 

The slope of our TF relation is in good agreement with that obtained from the calibration of the TF relation 
by G97 using SCI galaxies. However, we find a significantly larger TF scatter (cttf = 0.44) than the average scatter 
found by G97 {otf = 0.36), and a significantly larger zero point (Atf = —5.89 vs Atf = —6.09). We thus conclude 
that the scatter in the TF relation for field galaxies is larger and the zero point is smaller than the correspoding 
values for galaxies in clusters, confirming the results of D98 and Freudling et al. (1999). Our analysis also confirms 
that the TF scatter decreases with increasing galaxy rotation velocity, but this does not affect our estimates of /3/. 
We have also found that the velocity noise, <7„, increases with distance up to ~ 500 kms -1 and is significantly larger 
than the values obtained by VI and V2. Although the meaning of such a large a v is not clear, it is reassuring that it 
does not affect the general result of the VELMOD analysis and, in particular, the value of /3i. 

In summary, we have found that the PSCz velocity model provides a satisfactory fit to the velocities of SFI 
galaxies when a value of f3i = 0.42 ± 0.04 is used. This value is in good agreement with most of the recent v-v 
analyses that use the IRAS gravity field, inferred either from the 1.2Jy or the PSCz surveys, irrespective of the type 
of distance indicator used. Indeed, values of (3i in the range 0.4 - 0.6 are found using the TF relation for spiral galaxies 
(Davis et al. 1996, VI, V2, D98), the D n — a relation for early type galaxies (Nusser et al. 2000), Type la Supernovae 
(Riess et al. 1997) and the "surface brightness fluctuations" method in nearby galaxies (Blakeslee et al. 2000). Our 
estimate of /3j is also consistent with the results of v-v analyses which use a model gravity field derived from optical 
galaxies in the ORS survey (Santiago et al. 1995, 1996), once the stronger clustering of optical galaxies is taken into 
account. 

These values of f3j from v-v analyses are significantly smaller than those found using d-d methods (e.g. Sigad et 
al. 1998), even when both are applied to identical datasets. Further support for high values of (3 comes from power 
spectrum analyses of various velocity catalogues (Zaroubi et al. 1997, Freudling et al. 1999, Zaroubi et al. 2000), 
although in this case it is possible that accounting for nonlinear effects may help reduce the discrepancy (Silberman 
et al. 2000). Since v-v methods generally rely on the deviations of a set of observed velocities from a set of model 
velocities, differences between the two sets do not contribute to f3 but to the random errors, whereas this is not the 
case in power spectrum analyses. 

A different explanation for the discrepancy between v-v and d-d methods is that the v-v analyses are affected 
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by non-linear effects. All v-v methods implemented so far assume linear biasing and most of them also use linear 
gravitational instability theory. Recently new methods have been introduced to measure the nonlinearity in the bias 
relation in redshift surveys (e.g. Matarrese, Verde & Heavens 1997, Szapudi 1998, Sigad, Branchini & Dekel 2000, 
Feldman et al. 2000) and some of these have been applied to the PSCz survey. For example, Branchini et al. (in 
prep.) have applied the technique of Sigad et al. (2000) to the PSCz sample. Their preliminary results indicate that 
deviations from linear biasing are small and manifest themselves mainly as a anti-bias in low density regions. The 
amplitude of these effects can only lead to a modest increase of Pi, of the order of ~ 10%. Nonlinear motions might 
also change the value of /?/. However, Shaya, Peebles & Tully (1995) used a fully nonlinear model of the velocity field 
and still recovered a low value of /?. Similarly, VI performed a nonlinear VELMOD analysis in an attempt to break 
the degeneracy between fi m and bi and found indirect evidence that nonlinear motions are already accounted for in 
the model velocity field when this is predicted from a smooth density field. These results suggest that a treatment 
of nonlinear effects and the inclusion of nonlinear biasing prescriptions in v-v analyses will not change the inferred 
value of (5 significantly. Thus, it seems unlikely that the disagreement between the two methods will be eliminated 
purely by improving the existing v-v techniques. 
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Table 1. VELMOD analysis of mock catalogues. The top row lists the value of the true parameters. The second row lists the 
mean value, the error on the mean and, in parenthesis, the typical error in a single realization from VELMOD. Column 1: TF 
zero point ; column 2: TF slope; column 3: TF dispersion; column 4: ID velocity dispersion (in kms -1 ); column 5: value of 0i. 



Atf i>TF &tf &v Pi 

-6.1 7.33 0.36 150 1.0 

-6.096 ± 0.005(0.02) 7.31 ± 0.04(0.16) 0.36 ± 0.002(0.01) 150 ± 7(33) 1.01 ± 0.009(0.04) 



Table 2. VELMOD analysis of Mark III subcatalogucs. The corresponding results of the VI and V2 analyses are quoted in 
parenthesis. Column 1: Subcatalogue; column 2: TF zero point; column 3: TF slope; column 4: TF dispersion; column 5: ID 
velocity dispersion (in kms -1 ); column 6: f3j value and its 1-cr error 



Catalogue Atf ^tf <^tf <?v Pi 

A82 -5.91 (-5.96) 10.43 (10.36) 0.468 (0.464) 159 (125) 0.50 ± 0.07 (0.49 ± 0.08) 

MAT -5.84 (-5.75) 7.06 (7.12) 0.451 (0.453) 180 (125) 0.42 ± 0.10 (0.50 ± 0.11) 

MAT2 -5.83 (-5.80) 6.96 (7.16) 0.433 (0.430) 220 (130) 0.44 ± 0.07 (0.52 ± 0.05) 



Table 3. VELMOD analysis of various SFI subsamples. Column 1: rcdshift range (in kms -1 ) and additional cuts; column 2: 
number of galaxies column 3: TF zero point; column 4: TF slope; column 5: TF dispersion; column 6: ID velocity dispersion 
(in kms -1 ); column 7: p m in and its 1-cr error 



Sample 


N gal. 


Atf 


bTF 


<TTF 


0~ v 


Pi 


cz LG < 6000 


989 


-5.89 


7.19 


0.439 


250 


0.42 ± 0.04 


cz LG < 6000, r; > -0.1 


640 


-5.89 


6.74 


0.382 


230 


0.43 ± 0.04 


cz LG < 6000, r; < -0.1 


349 


-5.87 


7.04 


0.498 


248 


0.39 ± 0.04 


cz LG < 4000 


493 


-5.86 


7.32 


0.415 


216 


0.45 ± 0.05 


cz LG < 2000 


158 


-5.73 


6.67 


0.454 


206 


0.42 ± 0.09 


2000 < cz LG < 4000 


335 


-5.88 


7.41 


0.410 


196 


0.44 ± 0.07 


4000 < cz LG < 6000 


496 


-5.85 


7.54 


0.420 


535 


0.40 ± 0.09 
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